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Abstract  :  This  paper  presents  a  method  for  solving  a  stochastic  version  of  the  dynamic 
traffic  assignment  problem.  It  shows  that  a  globally  optimal  solution  may  be  obtained  by  a 
sequence  of  linear  optimizations.  A  decomposition  algorithm  for  this  procedure  is  presented 
that  efficiently  solves  large-scale  problems.  Solution  examples  with  up  to  sixty-six  thousand 

variables  are  described.  \ 
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1.  Introduction 


The  problem  of  determining  flows  in  a  network  to  optimize  an  objective  or  satisfy 
equilibrium  conditions  has  been  studied  in  a  variety  of  frameworks.  The  static  problem  has 
been  extensively  studied  (see,  for  example,  Dafermos  and  Sparrow  [3]  and  Potts  and  Oliver 
[7]),  but  the  dynamic  problem  has  received  less  attention.  Merchant  and  Nemhauser  [5] 
gave  a  mathematical  programming  formulation  of  this  dynamic  traffic  assignment  problem. 
They  showed  that  the  resulting  nonlinear,  nonconvex  program  has  a  piecewise  linear  version 
such  that  a  feasible,  globally  optimal  solution  exists  among  the  set  of  linear  program  optima. 

Ho  [4]  showed  that  a  global  optimum  could  be  obtained  by  solving  a  sequence  of  linear 
programs.  The  number  of  such  programs  is  at  most  N  +  I  where  N  is  the  number  of 
periods. 

The  previous  studies  have  all  assumed  that  the  exogenous  flows  into  nodes  are  known 
for  all  periods.  This  assumption  is  not  generally  true.  Instead,  the  flows  are  unknown 
and  represent  random  rates  of  arrivals  into  the  system.  The  stochastic  dynamic  traffic 
assignment  problem  is  then  to  determine  flows  along  the  arcs  to  minimize  a  convex  function 
of  traffic  congestion.  Flow  decisions  in  one  period  depend  on  the  realization  of  inputs  from 
previous  periods.  The  decisions  seek  to  minimize  expected  costs  over  future  periods. 

In  this  paper,  we  formulate  a  piecewise  linear  version  of  the  stochastic  dynamic  traffic 
assignment  problem.  This  representation  is  analogous  to  those  given  by  Ho  and  Merchant 
and  Nemhauser.  In  this  multistage  model,  uncertainties  are  resolved  period  by  period 
as  time  progresses.  We  show  that  a  successive  linearization  procedure,  similar  to  Ho’s, 
yields  a  globally  optimal  solution.  We  also  present  an  algorithm  based  on  the  multistage 
decomposition  method  of  Birge  [  1  j  to  implement  this  procedure.  Experimental  results  for 
this  algorithm  on  problems  with  up  to  sixty-six  thousand  variables  are  given. 


2.  Multistage  Stochastic  Program  Formulation 


The  multistage  stochastic  program  assumes  that  decisions  taken  in  period  i  depend 
on  outcomes  and  decision  made  in  periods  1, . . . —  1.  Period  «  decisions  cannot,  however, 
depend  explicitly  on  future  outcomes  (i.e.,  they  are  nonanticipative) .  In  every  period  i,  a 
finite  number  of  scenarios  represent  different  sets  of  inputs.  Each  scenario  l{  has  an  asso¬ 
ciated  probability  pi. .  For  every  scenario  in  period  t(t  >  0),  there  corresponds  a  scenario 
a(li)  in  period  »  -  1  which  is  the  immediate  ancestor  scenario  of  1,.  Each  /j  also  has  a  set 
A(li)  of  ancestor  scenarios  (one  in  each  period  *'  <  *)  that  is  defined  recursively  as  the  set  of 
currently  defined  ancestor  scenarios  and  their  immediate  ancestors.  There  also  exists  a  set 
of  immediate  descendant  scenarios  Di+l(li)  C  {1, ...,  Li+l},  in  period  «  +  1  for  every  t.  The 
descendant  scenarios  partition  {1,. . . ,  Li+1}  so  that  Ui,e{i  l{)  Pi+i(*i)  =  {1,-  ••  ,£,+  i} 
and  Di+ i (*)  n  A+i(i)  =  0  if  i  ^  j.  We  let  D(U)  be  the  set  of  descendant  scenarios  consist¬ 
ing  of  immediate  descendant  scenarios  of  »  and  all  future  (recursively  defined)  immediate 
descendants  up  to  the  planning  horizon. 

The  flow  decisions  are  i.yi  for  arc  j  in  period  »  and  scenario  l.  Other  notation  is 
consistent  with  Ho  with  decision  variables  X*}l  representing  weights  on  grid  points  c*  such 
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xi3  i  =  ^  cy  ■ 


For  period  «  =  0,  we  define  L0  =  1  and  let  x0]  and  A*;  be  interchangeable  with  i0;  i  and 
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G  =  (>/,£)  a  directed  graph; 


M  —  set  of  nodes  of  G ; 


S  =  set  of  arcs  (directed  edges)  of  G\ 
N  =  planning  horizon; 
i  =  index  of  time  period; 
j  =  index  of  arc  in  £  \j  =  1,  ...,o; 


g  =  index  of  node  of  M\q  =  l,...,n; 


n  =  index  of  destination  node; 

A(g)  =  {j  €  £  |  arc  j  leaves  node  q}\ 

B{q)  =  {j  €  £  |  arc  j  enters  node  ?}; 

Fu(q)  ~  external  input  at  node  q  under  scenario  /; 

*<» 

hufan)  =  co8t  of  xm  =  Y  hiixiit\ 

k  =  1 

diji  =  amount  of  traffic  admitted  to  arc  j  in  period  t  under  scenario  /; 

9i(x*i>)  =  amount  of  traffic  to  exit  from  arc  j  in  period  »  under  scenario  l 
KU ) 


We  make  the  following  assumptions  as  in  Ho[4]. 


(Al)  The  arcs  are  not  explicitly  capacitated. 


(A2)  Saturation  occurs  for  large  enough  x,  i.e.,  =  0. 

(A3)  The  nonnegative  slope  of  the  piecewise  linear  approximation  to  g}  is  strictly 


decreasing. 


(A4)  The  cost  functions  are  structured  so  that 


t. 


(a)  0  <  Afy  <  h*+ 1  for  all  i,j  and  k; 

(b)  If  is  closer  to  the  sink  than  j,  along  any  simple  path,  then  <  h\}i  for  all 


We  add  an  assumption  from  Merchant  and  Nemhauser  [5]: 

(A5)  0  <  ffj(x)  <  x. 

The  stochastic  dynamic  traffic  assignment  problem  (SDTAP)  is  then  formulated  as: 


SDTAP.  0 


■"in  £*=l +E2,lES,p.{Ef.T<i{,A{„}  =  o»*,W 

subject  to 

E™  c*a;„  =  EfJ’,1  W  W,  +  *>;.J  =  1 =  1 W,  SDTAP.  1 

Efii1  =  SSfW  -  »,*)*?»  +  <W  = 1 . *  - 1; 


i 

=  1... 

•  ,o;/'  €  Di+i(l);l 

—  1, . . . ,  Li ; 

^l€A(q) 

doy  = 

■  v*  <;) 

'  Z_ifc=i 

9*a  o,-,v<i#  n; 

^y€i4(fl) 

d»y,  = 

Fn(q) 

i  v-i 

f  <L>}'eB(q)  lsk=  1 

t 

=  i,.. 

i  c, 

A*,  = 

: 

=  1,  •••><*; 

r  do, 

>0,4, 

>0, 

» =  i, 

a,/  =  1, . . . , 

K} 

>0, 

j  =  i, 

, a,  Ar  =  1,...,. 

R’O);  \ 

K 

,  >0, 

»'  =  i, 

=  1,..., 

a,  / 

/=  i, 

. . . ,  Li,k  =  1, . . . , 

*0); 

1  xk 

Aoy 

=  1,12 

(;)-i 

=  i 

A,-,-,  = 

i;»  =  i 

=  =  1; 

/  ^oy 

>  0 

for  at 

most  two  consecutive  k, , 

1  Afyi 

>  0 

for  at 

most  two  consecutive  k.  ' 

SDTAP.  2 
SDTAP.  3 

SDTAPA 
SDTAP.  5 

SDTAP. 6 

SDTAP.l 

OSP 


SDTAP  is  a  multistage  stochastic  program  with  recourse  (see  Wets  [8]).  Decisions  A*, 
and  diji  are  made  knowing  P,>,  for i'  <  »  with  certainty  but  without  knowing  the  realizations 
of  future  inflows  P,.,  for  »'  >  ».  The  future  inflows  represent  the  discrete  realizations  of  the 
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random  external  input. 


The  ordered  solution  property  constraint  (OSP)  makes  SDTAP  a  nonconvex  program. 
Additional  optimization  must,  therefore,  be  performed  to  solve  it.  In  Ho  [4],  it  is  shown  that 
an  optimal  solution  can  be  found  for  the  deterministic  problem  by  successively  optimizing 
at  most  IV  +  1  objective  functions  for  the  linear  programming  part  of  SDTAP.  We  show 
that  a  similar  successive  optimization  approach  yields  a  globally  optimal  soultion  for  the 
stochastic  problem. 

The  difference  between  the  procedure  followed  below  and  Ho’s  procedure  is  our  order 
of  optimization  and  choice  of  objectives.  These  objectives  are  used  on  individual  scenarios 
to  maintain  the  form  of  the  stochastic  program.  We  first  state  a  result  from  Lemmae  A,B,C 
of  Ho  and  Lemma  1  of  Merchant  and  Nemhauser. 

Lemma  1.  Let  y  =  {A£,, <£>•<}  be  a  feasible  solution  of  SDTAP.1-7  that  violates  OSP  for 
»  =  r,  j  =  8,  and  l  —  t.  There  exists  a  feasible  solution,  y  =  {A*,,*/,,/},  to  SDTAP.1-7 
that  differs  from  y  only  for  j  =  a  when  t  =  r  and  l  =  t,  and  for  arcs  j  on  paths  from  s  to  n 
for  i  >  r  and  /  €  Pi(t).  The  solution  for  A  in  (1),  xr,t  =  xr,t,  the  solution  given  by  A.  The 
solution  y  satisfies  OSP  for  i  =  r  and  j  —  a  and: 

a. )  For  all  q  €  M,  i  =  0, . . . ,  N,  l  =  1 ,Lit  the  total  flow  reaching  q  on  or  before  t  is  at 

least  as  great  for  y  as  for  y. 

b. )  For  any  scenario  /  G  Ls  ,  the  scenario  l  objective  value, 

N  a  K(j) 

<4+i( a)  =  E  E  EE > g‘n+a a). 

»=0  l'eA(t)  ;  =  1  fc=  1 

c. )  The  solution  { \i]t , d,y( }  is  the  same  for  any  value  of  {A i'  >  *. 

Proof:  Applying  Lemma  1  of  Merchant  and  Nemhauser  on  each  scenario  /  =  1, . . . ,  Ln  and 


'.LI  L  I.u  wmwwpp 


its  ancestor  scenarios,  A{1),  directly  yields  the  results  above  except  for  the  nonanticipativity 
requirement  in  SDTAP.  We,  therefore,  need  to  obtain  the  same  A*  Y,  and  d,j/>  for  all  l  E  D(l'). 
By  the  construction  in  Merchant  and  Nemhauser,  (c)  above  holds  so  that  nonanticipativity 
is  maintained.! 

The  following  corollary  is  the  basis  for  the  successive  optimization  procedure. 


Corollary  1.  Consider  a  feasible  solution  y  =  {A*}l,di}t}  of  SDTAP.1-7.  Suppose  y*  = 
{A**, ,  d*  .| }  is  an  optimal  solution  of 


max  G,</<  (A) 

—  V'“  -k  \  fc 

—  1-1  i=\  Lj*=\  9}*i'ji' 

subject  to 

(SDTAP.  1  -  SDTAP. 7) 

(A 

=  (%>4 ■»);*  <  *v.Vj 

<G'N+l(A),/=l,...,L, 

(SUB  -  i'l') 


Then,  y*  has  no  OSP  violations  in  period  i'  and  scenario  l'. 

Proof:  Suppose  an  OSP  violation  occurs  in  y*  in  period  »' ,  scenario  /' ,  arc  f .  By  Lemma 
1,  there  exists  a  solution  y  =  {A,  d}  that  satisfies  OSP  on  arc  j‘  in  scenario  V  at  period  «' 
and  that  satisfies  the  constraints  in  (SUB-tT).  By  strict  concavity  in  Assumption  A3, 


K(j')  _  K(i') 

y.  9) 'K'j'i'  >  9j'^i')'e  ■ 

k=l 


(2) 


fc=i 


Hence,  G,>j<(A)  >  Gj-j«(A*),  a  contradictions  of  A*’s  optimality.! 

This  corollary  can  be  applied  repeated ly  from  period  1  to  N  for  each  scenario  to  obtain 
a  full  OSP  solution.  The  repeated  direct  solution  of  the  full  multistage  problem  using  stan¬ 
dard  linear  programming  procedures  may,  however,  be  too  computationally  burdensome 
to  implement.  The  structure  of  SDTAP  can,  however,  be  exploited.  In  [lj,  a  multistage 
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nested  Benders’  decomposition  method  (NDSP)  is  described  that  significantly  reduced  solu¬ 
tion  times  from  simplex  method  solutions  on  a  set  of  practical  test  problems.  NDSP  applies 
directly  to  SDTAP.1-7.  Obtaining  an  OSP  solution  using  objectives  Gw  requires,  however, 
the  addition  of  the  intertemporal  constraints,  G‘N  +  1(X)  <  GlN  +  i (A).  These  constraints  are 
not  directly  compatible  in  the  decompoation  procedure.  The  development  below  shows  how 
to  incorporate  these  constraints  into  a  procedure  to  obtain  OSP. 

To  simplify  notation  in  the  following  development,  let  y,(  =  {A*.,,  dij(},  let  SDTAP.l 
be 

N  Li 

Gjv+i(y)  =  (3) 

•  =0  1=1 

let  the  constraints,  SDTAP.l, 3,  and  5,  be 

A)S/o  =  &o,  (4) 

let  the  constraints,  SDTAP.2  and  4,  be 

£«-iV*-i,r  +  =  &,V/  e  =  1  =  1, . . . , L,_i;  (5) 

and  note  that  below  we  use  j  and  k  as  general  indices  that  are  not  restricted  to  denoting 
arcs  and  linear  segments. 

The  basic  (i,/)-subproblem  solved  in  NDSP,  given  some  solution,  y<_i,0(j)t  is  then 


„(,))  =  min  htyu 

+9n 

(NDS  -  1) 

subject  to  Ai  yu 

—  Gil  -  -B* — i  y» — i .  o  ( i ) ; 

[NDS  -  2) 

D«y<i 

—  ^il  >  r  ~  •  ■  •  >  ! 

( NDS  -  3) 

E;,yu 

+0« 

>  e'(, s  =  1,. . .  ,5,,; 

( NDS  -  4) 

Vu 

>  o. 

( NDS  -  5) 

Constraints  (NDS-3)  are  feasibility  cuts  to  maintain  a  feasible  solution  in  descendant  sce¬ 
narios  /.  The  (NDS-4)  constraints  are  optimality  cuts  on  6,t .  They  represent  an  outer 
linearization  of  Gn+i  restricted  to  descendants  of  /.  NDSP  repeatedly  solves  subproblems 
NDS  to  obtain  (yn,0u)  until  the  optimality  condition, 

PiAi  >  5Z  Pi+U'Qi+irO/ii),  (6) 

is  achieved  for  all  i,l.  (Note  that  pu  =  52i'eDi+l{t)  P»+ 

For  the  algorithm  below,  we  also  define 

N 

PuGN  +  i(y)  =  Y2  H  Wt'hi'Qet’,  (7) 

i'=n'eP(i ) 

as  the  contribution  to  the  objective  of  l  and  its  descendant  scenarios.  The  following  algo¬ 
rithm  obtains  OSP  in  all  t  and  /. 

Nested  Decomposition  of  Stochastic  IVaffic  Assignment  Problem  (NDSTAP) 

0.  Initialization.  Use  NDSP  to  obtain  y,  a  solution  to  SDTAP.1-7.  Let  »  =  —  1,  /  =  0.  Let 
Lq  =  1. 

1.  Check  for  OSP.  Let  1  =  1  +  1.  If  *  >  iV,  stop;  y  is  a  full  OSP  solution. 

If  /  >  Li,  let  t  =  i  +  1,  /  =  0,  return  to  1. 

Else,  check  OSP  in  y,j.  If  OSP  is  satisfied,  return  to  1. 

Else,  go  to  2. 

2.  Setup  for  OSP  at  (i,l).  Let  y<<(<  be  as  in  the  NDSP  solution  for  all  »'  =  i  +  1, . . . ,  N, 
and  /'  €  P(l).  For  each  such  subproblem  (»',/'),  let  the  NDS-3  and  NDS-4  constraints  be 
those  defined  at  the  end  of  Step  0.  Let  Kw  =  0  for  all  i'  =  i  + 1, . . . ,  N\ /'  €  P(l).  Go  to  3. 
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8.  Maximize  flow  at  i,l.  Solve  the  subproblem: 


max  Gu(yu) 

(8.0) 

subject  to  A+yu 

—  £il  —  1  S/i  —  l,a(f)  j 

(8.1) 

Kvh 

-z« 

=  -  Ey=o  -C('eX(i) 

(8.2) 

Kyu 

+#« 

-  GN+i(y)’ 

(8.3) 

Dri,y« 

>  dr,,r  =  1, . . . , Ru\ 

(8.4) 

Kyu 

+0il 

—  eit>  8  =  1»  4  •  •  ,  Siti 

(8.5) 

D-,y.i 

+$li  Zit 

>  d*2 ,k  =  1  ,...,Kui 

(8.6) 

y« 

>0; 

(8.7) 

to  obtain  (y,(, i,().  Let  i'  =  i  +  1.  Go  to  4. 

4-  Cheek  feasibility  at  i' .  If  i' 

>  N,  go  to  1. 

Else,  for  all  V  6  D{1)  at  solve 

min  hayin' 

+0i>i> 

(9.0) 

subject  to  Ai'yi'f 

=  & i>  -  .B, 

(9.1) 

hi-yi'i' 

—Zi'r 

(9.2) 

>  r  =  1 

(9.3) 

E-.^yi'i' 

^  ei'i'  >  ®  =  !>•••>  Sin' ; 

(9.4) 

Di'yyi'i' 

+d%'i'Zi<i' 

>  d*,2  ,k  =  1, . . . ,  Km' ,  i'  <  N: 

;  (9.5a) 

-Zi'f 

(9.55) 

Vi'V 

>o, 

(9.6) 

to  obtain  a  new 

If  all  subproblems  (9)  at 

i  are  feasible,  let  i'  =  »'  +  !,  return  to  4. 

Else,  for  (9)  infeasible  at  (*',/'),  add  a 

feasibility  cut,  (9.5a),  to  (*'  —  l,o 

(/')),  where 

(101) 
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(10.2) 


—  7(2 '< 

5*'  -!,<»(/')  =  » 3di'i>  +  n  4e,/(-  +  7T5 


(10.3) 


w  =  di*v  if  i'  <  N;  w  =  -G‘N+l( y)  if »'  =  N,  k  =  +  1,  and  =  1, ...  ,5,  are 

the  multipliers  corresponding  to  constraints  (9.j),  j  =  1, . . . ,  5,  respectively,  that  obtained 
the  infeasibility  condition.  Update  for  each  infeasible  (s', /')  subproblem  (9).  Let 

t'  =  »'  —  1.  If  »'  =  »,  go  to  3.  Else,  return  to  4. 


The  NDSTAP  steps  1  to  4  differ  from  NDSP  iterations  in  the  use  of  the  linking  con¬ 
straints  (8.2)  and  (9.2).  NDSTAP  also  does  not  add  any  optimality  cuts  (8.5)  and  (9.4) 
because  (9.5b)  guarantees  optimality. 

Other  strategies  can  also  be  used  in  proceeding  between  periods.  The  above  strategy 
solves  all  problems  at  one  period  proceeding  to  the  next  period.  An  alternative  is  to  solve 
all  descendant  subproblems  for  the  current  scenario  at  »'  before  proceeding  to  the  next  i' 
scenario.  The  above  implementation  required  few  iterations,  however,  and  was  used  for  the 
computational  tests.  Convergence  of  the  NDSTAP  iterations  to  an  OSP  solution  is  given 
in  the  following  theorem. 


Theorem  1.  NDSTAP  terminates  in  a  finite  number  of  iterations  of  Steps  1  to  4  with  a 
solution  {yit}  that  satisfies  OSP  for  all  i  and  /. 


Proof:  OSP  is  obtained  for  i  =  0  if  Steps  2  to  4  terminate  for  »  =  0,  according  to  Corollary 
1.  Given  OSP  at  t  =  0,  Corollary  1  and  induction  establish  OSP  for  all  i  and  /  if  Steps  2 
to  4  terminate  finitely  for  each  non-OSP  solution  found. 

To  show  that  Steps  2  to  4  terminate  finitely,  note  that,  because  the  constraints  (8.3- 
5)  and  (9.3-4)  are  outer  linearizations,  any  feasible  solution  {t/<;}  such  that  G;v  +  i(y)  < 


r.Mj 


On +i{y)  must  satisfy  (8.3-5)  and  (9.3-4)  for  all  t  and  /.  These  constraints,  therefore,  do 
not  exclude  any  solutions  y  that  satisfy 

G‘N+1(y)  <  0‘N+l(y)  for  all  /  =  1, . . . ,  Ln  ■  (11) 

Since  (8.2)  and  (9.2)  are  always  satisfiable,  we  must  establish  that 

1.  there  exists  yit  and  |  t'  >  t,l'  £  0(1)}  that  are  feasible  for  (8.1-5)  and  (9. 1-4, 5b) 
for  all  »', 

2.  Constraints  (8.6)  and  (9.5a)  do  not  exclude  any  feasible  y  that  satisfies  (11); 

3.  Only  a  finite  number  of  constraints  (8.6)  and  (9.7)  are  generated. 

The  first  condition  is  guaranteed  by  Corollary  1  as  stated  above.  Condition  2  is  true  by 
induction.  Assume  that  v  =  yjv-i, <,(/'),  u  =  ZN-i,a(i')  is  input  to  subproblem  (9)  for  ( N,l '). 
Then,  there  exist  multipliers  x1 , . . . ,  x5  (x3  =  x4  =  0  since  (9.3)  and  (9.4)  are  vacuous  for 
i  =  N),  such  that 

x1Ajv  -f-jr2hjv  <  0, 

-x2  -x5  =  0, 

x5  >  0,  (12) 

and 

x1^ vv  ~  xxBS-\v  -  x2u  -  x5Gltf  +  l(y)  >  0.  (13) 

In  order  for  t>  and  u  to  be  feasible,  therefore,  we  require 

*1Bn-iV  +  x2u  >  x'Zni1  ~  xbGlN  +  l(y).  (14) 

Constraint  14  is  the  feasibility  cut  (9.5a)  as  defined  in  (10).  Hence,  no  constraints  (9.5a) 
for  i'  =  N  -  1  exclude  any  feasible  y.  By  induction,  the  same  argument  applies  for  all 
constraints  (8.6)  and  (9.5a),  establishing  Condition  2. 


Condition  3  ia  true  by  induction  also.  For  i  =  N,  only  extreme  solutions  of  (12)  need 
be  considered.  Therefore,  only  a  finite  number  of  cuts  (14)  are  generated.  Hence,  only  a 
finite  number  of  constraints  (9.5)  are  added  to  subproblems  for  i'  =  JV  —  1.  By  induction, 
this  is  true  for  all  »'.■ 

3.  Computational  Results 

The  algorithm  NDSTAP  was  coded  in  FORTRAN  as  an  extension  of  the  NDSP  code 
([l]).  The  code  applies  to  problems  with  up  to  three  periods  ( N  =  2)  with  multiple  de¬ 
scendants  in  each  period.  Scenarios  with  single  immediate  descendants  can  be  grouped 
together,  however,  to  solve  problems  with  N  >  2  and  L2  =  •  •  •  =  LN .  This  situation  corre¬ 
sponds,  in  practice,  to  the  ability  to  predict  future  inflows  with  certainty  given  observation 
of  the  first  three  periods  of  inflows. 

The  NDSTAP  code  takes  advantage  of  single-descendant  scenarios  in  the  last  period 
by  combining  periods  2  to  N  into  a  single  subproblem: 


min 

GjV-fc(y2Ji  •  •  •  i! Ini) 

(15.0) 

subject  to 

A2y2i 

—  -  #iyi,a(i); 

(15.1) 

BjiVjt  +Ay+Uy,+1,, 

=  Z,i,j  =  2,...,iV  — 1; 

(15.2) 

EJU  h>'V,i 

-  ~  ^iyi,a(J)  -  h0y o,a(a(/)); 

(15.3) 

Gn->(v 2i,---,ysi) 

>  GN-,,8  =  0,  ...,k  -  1(  if  k  >  0); 

(15.4) 

y a 

>0,j  =  2,...,AT; 

(15.5) 

where  Gy- 

»(y 3«  .  •  •  •  .  VNl)  =  Hf=2*  9}V)I  is 

the  total  flow  from  each  arc  in  periods  2 

to  N-s, 

and  GN is  the  value  of  Gy-,  obtained  for  subproblem  (15)  when  solved  with  objective 
Gy~,.  Subproblem  (15)  is  solved  instead  of  (8)  for  increasing  values  of  k  (starting  at  k  =  0) 
until  an  OSP  solution  is  obtained  for  (y2 i,-,yyi)-  This  is  the  successive  optimization 
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procedure  in  Ho,  which  again  achieves  OSP  finitely.  Its  advantage  over  solving  (8)  and 
(9)  repeatedly  is  that  a  single  solution  of  (15)  may  obtain  OSP  for  all  y,<,«  =  2, . . . ,  N, 
whereas  several  subproblems  (8)  may  be  required.  The  use  of  a  similar  objective  with  several 
descendant  scenarios  (i.e.,  when  Li  <  Li+1),  however,  still  requires  several  optimizations 
and  is  not  advantageous  over  (8)  and  (9).  The  strategy  in  NDSTAP  requires  only  one 
additional  optimization  to  obtain  OSP  in  period  0,  which  may  be  the  only  decision  actually 
implemented  before  another  problem  is  solved. 

NDSTAP  takes  further  advantage  of  the  structure  of  (15)  by  checking  for  an  OSP  bas  is, 
i.e.,  a  basis  of  (15)  such  that  A* ,  is  basic  for  at  most  two  consecutive  k.  When  an  OSP 
basis  is  found,  we  need  only  check  its  feasibility  in  other  subproblems  (15)  for  different  l  to 
find  an  OSP  solution  for  these  scenarios.  In  practice,  this  quickly  produces  OSP  solutions 
for  all  scenarios. 

The  test  problems  were  generated  with  the  same  structure  as  Ho’s  test  problems  in 
[4],  using  Ho’s  TAPGEN  generator  to  create  a  deterministic  problem  and  the  STGEN 
generator  used  in  [l]  to  create  the  stochastic  problem.  All  examples  use  the  network  in  [4] 
with  seven  nodes  and  twelve  arcs.  They  also  have  five  periods  (N  =  4). 

The  test  problems  include  varying  numbers  of  scenarios  and  different  period  0  input 
flows.  The  scenario  numbers  (Lj  =  L3  =  LA)  range  from  six  to  eighty-one.  The  period 
0  inputs  are  characterized  as  “high,”  “medium”  or  “low.”  Each  test  problem  has  three 
period-one  scenarios  {Li  =  3).  The  input  flows  are  summarized  in  Table  1.  The  period- 
two  scenarios  are  constructed  from  combinations  of  the  period-two  values  in  Table  1  for 


Table  1.  Input  Flows 


NODE 

PERIOD 

0 

1 

Input  Level 

Scenario 

High 

Medium 

Low 

1 

2 

3 

1 

10 

30 

50 

10 

15 

30 

2 

10 

30 

50 

20 

30 

30 

3 

10 

30 

50 

5 

10 

30 

4 

10 

30 

50 

0 

30 

30 

5 

10 

30 

50 

0 

30 

30 

6 

10 

30 

50 

0 

30 

30 

NODE 

PERIOD 

2 

3 

4 

1 

30 

20 

5 

2 

70 

40 

30 

3 

20 

10 

5 

4 

(10,20,30) 

30 

40 

5 

(15,25,35) 

60 

50 

6 

(0,15,30) 

10 

20 

The  size  characteristics  of  the  linear  programs  (SDTAP.1-7)  for  each  test  problem  are 
given  in  Table  2.  That  table  gives  the  number  of  period-two  scenarios,  the  number  of 
constraints  and  the  number  of  variables  for  each  SDTAP  solved. 


Table  2.  Test  Problem  Characteristics 


Problem  Numbers 

Scenarios 

Constraints 

Variables 

1-3 

6 

720 

5,952 

4-6 

12 

1,296 

10,800 

7-9 

24 

2,448 

20,496 

10-12 

81 

7,920 

66,552 

The  test  problems  were  solved  on  the  IBM  3033/4381  Computer  Network  at  the  Naval 


Postgraduate  School,  Monterey,  California,  using  the  VS  FORTRAN  compiler  under  the 
VM/CMS  operating  system.  The  results  from  solving  each  problem  are  given  in  Table 
3.  The  number  of  NDSP  “major  iterations”  refers  to  the  number  of  subproblems,  (NDS), 
solved  (excluding  multiple  solutions  of  last  period  subproblems  with  the  same  ancestor 
input).  NDSTAP  “major  iterations”  refer  to  the  number  of  subproblems  (8),  (9)  and 
(15)  solved.  “Simplex  iterations”  refer  to  the  number  of  simplex  iterations  performed  on 
subproblems  NDS,  (8),  (9)  and  (15).  The  CPU  second  times  (CPUs)  do  not  include  initial 
input  and  final  output  times  but  do  include  time  for  iteration  logging. 

The  number  of  OSP  violations  found  for  each  test  problem  in  Step  1  of  NDSTAP  are 
given  in  Table  4  by  period  (with  periods  2  to  4  grouped  together).  Whenever  an  OSP 
feasible  basis  in  (15)  was  found  for  some  scenario  /,  the  corresponding  OSP  solution  was 
used  for  that  scenario  so  that  no  OSP  violations  were  recorded  for  /. 


Table  3.  Test  Problem  Iterations  and  Times 


Problem 


NDSP 


1 

2 

3 


4 

5 

6 


7 

8 
9 


0 

1 

2 


Problem 


Period  0 
Input 


High 

Medium 

Low 


High 

Medium 

Low 


High 

Medium 

Low 


High 

Medium 

Low 


Period  0 
Input 


High 

Medium 

Low 


High 

Medium 

Low 


High 

Medium 

Low 


High 

Medium 

Low 


Major 

Iterations 

Simplex 

Iterations 

46 

326 

62 

346 

42 

353 

46 

402 

62 

449 

40 

364 

46 

468 

62 

535 

40 

480 

52 

656 

68 

1,024 

50 

729 

Major 

Iterations 


NDSTAP 


Simplex 

Iterations 


0.066 

0.059 

0.075 


0.103 

0.128 

0.162 


.140 

.188 

.290 


Table  4.  OSP  Violations 


I 


Problem 

Period 

Total 

0 

1 

2,3,4 

1 

2 

15 

19 

36 

2 

1 

12 

16 

29 

3 

0 

6 

20 

26 

4 

1 

13 

22 

36 

5 

1 

15 

18 

34 

6 

0 

7 

18 

25 

7 

1 

■ 

20 

35 

8 

1 

23 

39 

9 

0 

31 

38 

10 

1 

12 

19 

32 

11 

1 

16 

25 

42 

12 

0 

4 

34 

38 

Note  that  higher  input  flows  generally  produce  more  OSP  violations  (Table  4).  Ho 
also  observed  this  in  deterministic  problems.  Observe  in  Table  3,  however,  that  higher 
input  flow  problems  are  generally  easier  to  solve.  This  may  be  due  to  a  reduced  number  of 
feasible  bases  for  higher  input  flows  leading  to  fewer  iterations. 

The  results  in  Tables  3  and  4  indicate  that  NDSTAP  can  efficiently  produce  OSP 
solutions  for  large-scale  linear  programs  SDTAP.  The  test  problems  were  not  solved  by  a 
direct  simplex  method  implementation  because  time  and  memory  requirements  exceeded 
the  VM/CMS  limitations.*  The  true  advantage  of  NDSTAP  is,  indeed,  for  problems  that 

*  To  give  some  comparison,  however,  a  smaller  SDTAP  problem  (without  OSP)  with 
the  same  network,  511  constraints,  and  1318  variables  was  solved  in  [1],  NDSP  solved  this 
problem  in  0.76  CPU  seconds  on  an  Amdahl  V8  with  35  major  iterations  and  92  simplex 
iterations.  MINOS  (Murtagh  and  Saunders  [6])  solved  the  same  problem  on  the  same 
machine  in  10.37  CPU  seconds  with  337  simplex  iterations. 


cannot  be  solved  by  a  direct  simplex  method  implementation.  We  demonstrate  here  that 
OSP  solutions  can  efficiently  be  obtained  for  examples  of  these  difficult  problems. 

4.  Conclusion 

This  paper  presented  a  formulation  of  a  stochastic  dynamic  traffic  assignment  problem 
with  uncertain  input  flows  and  dynamic  decisions  given  knowledge  of  previous  inputs.  A 
successive  optimization  procedure  to  obtain  a  globally  optimal  solution  was  then  developed 
and  shown  to  converge.  The  efficiency  of  this  nested  decomposition  procedure,  NDSTAP, 
was  demonstrated  on  a  set  of  test  problems  with  up  to  almost  eight  thousand  linear  pro¬ 
gramming  constraints  and  over  sixty-six  thousand  variables.  An  NDSTAP  implementa¬ 
tion  solved  all  problems  in  under  five  CPU  seconds  on  an  IBM  3033/4381.  These  results 
demonstrate  the  practicality  of  decompostion  approaches  for  stochastic  traffic  assignment 
problems  and  other  general  multistage  stochastic  linear  programs. 
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